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SUMMARY 

Data  on  fracture  toughness  for  D6ac  steel  is  analysed  for  conformity  with  three 
probability  distributions.  The  three  parameter  extreme  value  distribution  is  selected,  and 
parameters  are  estimated  for  specified  conditions. 
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I.  INTRODUCTION 

'~-^A.R.L.  have  recently  applied  reliability  theory  to  structural  fatigue,  allowing  estimation 
of  safe  inspection  intervals  in  service  tflefs.  I awd'f).  The  method  recognises  various  inputs 
which  influence  the  risk  of  failure.  One  input  is  the  distribution  of  residual  strength  of  structures 
all  cracked  to  the  same  extent.  This  is  important  in  view  of  the  high  variability  in  fracture 
toughness  for  some  materials,  and  the  corresponding  variability  in  residual  strength  for  a given 
crack  size. 

The  feasibility  of  applying  the  method  was  investigated  using  the  example  of  an  aircraft 
having  many  components  made  from  D6ac  steel  heat  treated  to  strength  levels  of  220-240  k.s.i. 
Five  locations  in  D6ac  components  were  considered  possible  sites  for  fatigue  failures.  Much 
fracture  toughness  data  for  D6ac  steel  was  available.  This  was  analysed  to  identify  its  distri- 
bution, and  to  estimate  its  parameters  at  each  failure  site. 

This  report  presents  this  analysis. 


2.  GENERAL  REQUIREMENTS  ' 

The  requirement  is  to  identify  and  estimate  parameters  of  fracture  toughness  distributions 
for  the  five  failure  sites  at  temperatures  which  are  appropriate  for  two  types  of  condition.  These 
conditions  are  firstly  flight  operating  conditions,  and  secondly  low  temperatures  corresponding  to 
a cold  proof  load  test,  to  which  the  aircraft  is  subjected  before  entering  service.  These  tempera- 
tures are  specified  in  detail  in  section  4. 

It  should  be  noted  that  although  all  the  failure  sites  are  in  D6ac  components,  they  were 
not  all  heat  treated  by  the  same  process,  and  thus  mean  fracture  toughness  at  the  same  temperature 
could  be  expected  to  vary  among  the  failure  sites. 


3.  DATA 

Samples  of  fracture  toughness  data  are  given  in  Refs.  3 and  4.  Ref.  3 gives  extensive  data 
obtained  during  a multi-laboratory  test  programme  in  the  U.S.  Ref.  4 gives  a smaller  quantity 
of  dau  from  D6ac  components  from  the  aircraft. 

The  data  from  Ref.  3 covered  a wide  range  of  testing  parameters,  and  the  following  para- 
meters were  shown  to  have  an  important  effect  on  the  mean  and/or  variability  of  fracture 
toughneu  Kie'. 

(a)  Testing  Temperature. 

(b)  Local  quench  rate,  which  is  determined  by  the  heat  treatment  process  and  also  by 
parameter  (r). 

(r)  Thickness  of  material  during  heat  treatment.  (As  well  as  being  variable  itself  the  thickness 
during  heat  treatment  often  differed  from  the  thickness  of  the  specimen  tested). 

(<f)  Type  of  specimen.  Three  types  were  used;  compact  tension,  surface  flaw  and  double 
cantilever  beam. 

(e)  Materia)  batch  number. 

(/)  Position  of  specimen  within  the  part. 

Ref.  S gives  estinuites  of  mean  Xu  values  for  each  failure  site,  and  also  the  variation  of 
Xie  with  temperature. 

4.  ANALYSIS  PROCEDURE 

The  analysis  consisted  of  two  stages.  Firstly  it  was  assumed  that  a common  form  of  distri- 
bution for  A/,  of  06ac  steel  exists  under  all  testing  conditions,  and  this  form  was  estimated 
by  means  of  standardising  and  pooling  data  from  different  samples  to  give  the  largest  possible 
pooled  data  sample.  Secondly  the  parameters  of  the  distributions  were  estimated  by  considering 
data  relevant  to  etKh  failure  site  at  appropriate  temperatures. 
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4.1  Form  of  Distributkm 

Refs.  3 and  4 contain  a small  proportion  of  data  stated  to  he  invalid,  and  these  are  excluded 
from  the  present  analysis.  The  valid  data  were  grouped  into  10  samples,  each  having  a charac- 
teristic testing  temperature,  heat  treatment  and  specimen  type.  This  grouping  gave  samples 
which  included  variability  associated  with  batch  number,  geometric  variations  and  inherent 
variability  caused  by  micro  defects.  A total  of  406  data  points  were  obtained.  The  deuils  of 
the  samples  are  listed  in  Table  I.  It  is  noted  that  data  at  -40°F  and  -65°F  are  combined  into 
the  same  sample.  This  was  desirable  because  of  the  small  sample  sizes  at  each  temperature 
alone,  and  was  considered  justified  by  the  small  differences  in  data  between  the  two  temperatures. 

For  pooling  purposes  the  data  were  standardised  as 

x = (K  - R)IS 

where  R.  S represent  various  means  and  standard  deviations  for  the  relevant  subsamples.  A 
similar  variate  based  on  log  R was  used  to  compare  with  the  log  normal  distribution.  The 
standardised  data  were  then  pooled  into  a common  ordering  system.  The  empirical  distribution 
of  the  data  was  plotted  on  normal  probability  paper  (Figs.  I,  2)  to  compare  it  with  the  normal 
distribution.  It  should  be  noted  that  in  Figures  I and  2 every  point  has  been  plotted  at  prob- 
abilities greater  than  90%  and  less  than  10%,  whereas  between  these  values  every  tenth  point 
has  been  plotted.  An  inspection  of  the  common  ordering  system  showed  that  there  was  no 
strong  tendency  towards  clustering  of  data  from  any  particular  sample.  This  supports  the  initial 
assumption  that  a common  distribution  form  exists  for  all  samples,  and  that  data  from  various 
samples  can  be  pooled  after  standardising.  Figs  I and  2 indicate  that  the  Normal  and  log  Normal 
distributions  are  fairly  good  fits  to  the  data. 

The  data  from  group  5 (Table  I).  consisting  of  the  largest  sample  of  154  data  points,  was 
then  compared  with  the  two  parameter  extreme  value  distribution,  defined  as 

P{K)  - I - exp{ -(A/u)*} 
where 

K fracture  toughness  Ri^, 

P(R)  - probability  that  R is  less  than  the  nominated  value, 
u = characteristic  value  of  R,  for  which  P(R)  = 0-63, 
a - dispersion  and  shape  parameter. 

This  sample  alone  was  used  because  there  was  no  readily  available  method  of  standardising 
the  data  from  different  samples  which  is  compatible  with  the  extreme  value  distribution.  In  this 
instance  the  use  of  one  sample  alone  is  justified  by  its  large  size,  and  also  by  the  point  made 
in  the  previous  paragraph  that  there  was  no  strong  tendency  for  any  sample  to  appear  clustered 
in  the  total  pool  of  data.  Fig.  3 compares  the  data  with  the  two  parameter  extreme  value  dis- 
tribution, and  clearly  the  fit. 

As  an  alternative  the  three  parameter  extreme  value  distribution  was  investigated.  This  is 
defined  as 

where 

I = least  value  of  R,  for  which  P(R)  0.  and  v,  a are  as  before. 

The  result  is  shown  in  Fig.  4.  Clearly  this  is  an  improvement  on  the  two  parameter  distribution. 

Comparing  Figs  I,  2 and  4 suggests  that  the  Normal,  log  Normal  and  three  parameter 
exiromc  value  distributions  all  give  a fairly  reasonable  fit  to  the  data.  The  Kolmogarov-Smimoff 
test  (Ref.  6),  a non-parametric  test  which  makes  no  specific  assumptions  about  the  form  or 
I arameters  of  the  distribution  being  tested,  also  indicate  that  both  distributions  fitted  the  data 
adequately.  The  tests  were  significant  at  approximately  10%,  15%  and  20%  respectively  for 
the  Normal.  log-Normal  and  the  three-parameter  extreme  value  distributions.  Thus  it  cannot 
be  concluded  statistically  which  of  these  three  is  preferable. 

The  extreme  value  distribution  was  selected  because  it  allows  the  choice  of  a least  possible 
value  of  X/,;  the  normal  distribution  predicts  a finite  probability  of  achieving  a very  low  (or 
even  negative)  value  of  R/,.  In  practice  it  is  reasonable  to  expect  a non-zero  least  value  to  be 
imposed  by  manufacturing  quality  control. 
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4.2  Pwamten  of  DbtribvHoat 

F.stimates  of  the  parameters  u,  t and  « are  required  for  conditions  of  heat  treatment 
temperature  and  specimen  type  which  are  appropriate  for  the  five  failure  sites.  As  regards 
specimen  type,  the  data  in  Ref.  3 suggests  that  the  surface  flaw  type  gave  higher  values  of  fracture 
toughness  than  the  other  specimen  types.  For  this  reason  surface  flaw  specimen  data  were  not 
considered  when  estimating  parameters  of  distributions.  The  heat  treatment  processes  for  the 
various  failure  sites  were  known  in  terms  of  the  processes  used  for  the  data  in  Refs.  3 ai  d 4, 
As  discussed  above  in  section  2,  temperatures  were  required  to  be  appropriate  for  cold  proof 
load  test  conditions  as  well  as  for  normal  operating  conditions.  It  was  anticipated  that  tiie  proof 
load  test  could  be  applied  at  either  -40"F  or  -65‘'F,  and  both  these  temperatures  were  con- 
sidered Operating  temperatures  could  typically  vary  from  0°F  to  75"F,  with  short  periods  of 
high  speed  flight  at  higher  temperatures.  In  view  of  the  reduction  in  A;,  at  low  temperatures,  it 
was  decided  to  consider  operating  temperatures  of  0°F  and  20°F  as  well  as  75°F.  To  estimate 
the  distribution  parameters,  appropriate  samples  from  Refs.  3 and  4 were  used,  together  with 
dau  from  Ref.  5 giving  mean  Ki,  for  each  failure  site  and  the  variation  of  AT/,  with  temperature. 
It  is  noted  that  the  mean  values  in  Ref.  5 in  most  cases  differed  from  the  corresponding  sample 
means  from  Refs.  3 and  4 because  of  variations  in  local  quench  rate  within  the  part,  which  are 
to  some  extent  deterministic. 

The  procedure  was  to  select  data  samples  from  Refs.  3 and  4 having  the  temperatures 
(nearest  to  the  required  temperatures),  heat  treatment  and  specimen  type  for  each  failure  site. 
In  many  cases  these  distributions  were  shifted  to  give  the  appropriate  mean  values  from  Ref.  S, 
making  due  allowance  for  temperature  variations  (also  from  Ref.  5),  without  altering  the  shape 
of  the  distributions.  In  carrying  out  this  mean  shift  it  was  assumed  that  either  the  standard 
deviation  or  the  coefficient  of  variation  remained  constant.  Coefficient  of  variation  was  held 
constant  if  the  mean  shift  was  required  to  account  for  within  the  part  variations  of  AT/,  at  the 
same  temperature,  and  standard  deviation  was  held  constant  for  mean  shift  due  to  temperature. 
This  was  because,  for  the  data  in  Ref.  3,  the  mean  and  standard  deviation  were  both  reduced 
by  lowering  temperature,  whereas  the  standard  deviation  did  not  appear  to  vary  systematically 
with  sample  mean  at  the  same  temperature.  Finally  the  extreme  value  parameters  were  estimated 
graphically  using  extreme  value  probability  paper  and  the  expressions  in  Ref.  7 for  relation 
ships  between  the  various  parameters.  The  results  are  shown  in  Table  II. 

5.  CONCLUSIONS 

1 . The  analysis  suggests  that  a common  distribution  form  exists  for  the  fracture  toughness 
data  considered,  although  the  parameters  vary  greatly  according  to  the  testing  conditions. 

2.  The  daU  is  fitted  adequately  by  Normal,  log-Normal  and  three-parameter  extreme  value 
distributions.  This  generally  agrees  with  the  conclusions  of  reference  8 for  airframe  materials. 
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TABLE  I 

Dncrlpdoa  of  A/,  Data  SaaiH** 


(1) 

Sample 

Number 

(2) 

Specimen 

Type 

(3) 

Heat 

Treatment 

(4) 

Temp. 

(“F) 

(5) 

Sample 

Size 

(6) 

Source 
of  Data 

(7) 

Sample 

Mean 

(8) 

Sample 

S.D. 

1 

CT 

E 

-40* -65 

16 

Ref.  3 

36-7 

3-47 

2 

CT 

A 

70 

51 

Ref.  3 

95-8 

4-93 

3 

CT 

B 

70 

13 

Ref.  3 

85-6 

9-89 

4 

CT 

B 

-40  A -65 

13 

Ref.  3 

45*6 

3-34 

5 

CT 

E 

70 

154 

Ref.  3 

64-7 

11-20 

6 

DCB 

E 

70 

21 

Ref.  3 

65-9 

11-00 

7 

CT 

H 

70 

33 

Ref.  3 

56-5 

J3-64 

8 

SF 

E 

70 

57 

Ref  3 

81-2 

18-50 

9 

SF 

E 

Unknown 

20 

Ref  3 

78-7 

15-36 

10 

CT 

E 

70 

28 

Ref  4 

52- 1 

5-65 

(2)  — Compact  Tension,  Double  Cantilever  Beam,  or  Surface  Flaw. 

(3)  — "A"  and  “B”  — High  Toughness  Process 

“E”  — Medium  Toughness  process 
“H"  — Low  Toughness  Process 

(4)  — Dau  at  — 40® F and  — 65°F  was  combined. 
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FIG.  2 COMPARISON  OF  POOLED  DATA  WITH  LOG  NORMAL  DISTRIBUTION 
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FIG.  4 COMPARISON  OF  DATA  WITH  THREE  PARAMETER  EXTREME  VALUE  DISTRIBUTION 
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